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It has been suggested that codimension-two braneworlds might naturally explain the vanishing 

of the 4D effective cosmological constant, due to the automatic relation between the deficit angle 

and the brane tension. To investigate whether this cancellation happens dynamically, and within 

the context of a realistic cosmology, we study a codimension-two braneworld with spherical extra 

dimensions compactified by magnetic flux. Assuming Einstein gravity, we show that when the brane 

contains matter with an arbitrary equation of state, the 4D metric components are not regular at the 

brane, unless the brane has nonzero thickness. We construct explicit 6D solutions with thick branes, 

treating the brane matter as a perturbation, and find that the universe expands consistently with 

standard Friedmann-Robertson-Walker (FRW) cosmology. The relation between the brane tension 

^+ ' and the bulk deficit angle becomes A = 2-KGe{p — 3p) for a general equation of state. However, this 

C^ ' relation does not imply a self-tuning of the effective 4D cosmological constant to zero; perturbations 

c7^ , of the brane tension in a static solution lead to deSitter or anti-deSitter braneworlds. Our results 

04 ' thus confirm other recent work showing that codimension-two braneworlds in nonsupersymmetric 

Einstein gravity do not lead to a dynamical relaxation of the cosmological constant, but they leave 

O , open the possibility that supersymmetric versions can be compatible with self-tuning. 
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I. INTRODUCTION 



The impressive correspondence between the ever-increasing amount of high-precision cosmological data 0, Q and 

^O ' the theoretical predictions of the A-Cold Dark Matter (ACDM) model is undeniably one of the greatest success 

^D \ stories of modern physics. Casting a shadow over this success however is the nagging fact that there is to date no 

'^ • satisfactory understanding of just what the dark energy (A) component is, and why its properties are what they are. 

^r Working under the assumption that this term is a cosmological constant due to vacuum energy, naive estimates of 

,J^ . its size come out a staggering 60 — 120 orders of magnitude too large (the number depending on one's beliefs about 

I ' supersymmetry), which until recently led most physicists to suspect that there must exist some as-yet undiscovered 

Oh, mechanism which exactly cancels its effect. The fact that observations now tell us that this term is not exactly zero, 

^ ' but nevertheless extraordinarily small by particle physics standards, means that finding a theoretical framework where 

'^"1 ' we can understand its scale and origin has become a pressing challenge for the theoretical physics community. 

^ One promising line of attack on the cosmological constant problem, which has recently garnered renewed interest, 

centers around ideas which arose some twenty years ago with the work of Rubakov and Shaposhnikov '^. These 

authors realized that in certain six-dimensional models, the expansion rate of the three large spatial dimensions is an 

d ' integration constant and thus independent of the vacuum energy. Several papers followed, e.g. [J- Hi which studied 

features such as stability and cosmology in these models. 

Lately, the advent of braneworld scenarios _8|- JJJ has rekindled interest in these ideas as possible ways of addressing 
the cosmological constant problem [i3|-[l3- I^ the new perspective, the old six dimensional models can be recast 
as codimension-two braneworlds that have the interesting property that the brane tension, rather than leading to 
expansion, induces a deficit angle in the bulk. Thus the Hubble expansion measured by observers on the brane would 
be insensitive to the vacuum energy associated with any field theory confined to that brane. This has led to the hope 
that there might exist a mechanism whereby if the brane tension suddenly changed, as in a phase transition, the 
deficit angle would dynamically relax to a value which cancels the tension's contribution to the expansion rate, thus 
providing at least a partial solution to the cosmological constant problem. (It would be only a partial solution since 
it is still necessary to do one fine-tuning of bulk contributions to the effective 4D cosmlogical constant.) 

However there are open questions concerning this scenario. It is not at all clear that such a tuning mechanism is 
really possible within Einstein gravity [l3-|23- Furthermore, it has been argued that it is impossible to find solutions 
with any matter sources on the brane except for those with the equation of state of pure tension |l5l l2ll | , which is a 
major obstruction to building more general models that would allow one to recover FRW cosmology at late times. 
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Much effort has been devoted to studyin g th e cosmology of codimension-one branes, and the associated corrections 
to standard cosmology are well-known |22j-|24|. Lately codimension-two braneworlds have received increasingly more 
attention, but most of this work has focused on static solutions or issues of stability |23-|23' while the papers 
addressing cosmology have considered almost exclusively only pure tension branes (see, however, |23,|4^-|43 for 
recent proposals for generalizing the brane equation of state, and |44| for a different perspective on the matter) . 

In this paper, we argue that the apparent obstruction to finding solutions with a more general equation of state 
for codimension-two branes is not a problem in principle; rather it stems from an unrealistic expectation that the 
bulk geometry should be regular in the vicinity of the brane. This singularity happens to be absent for branes whose 
equation of state is that of pure tension, p — —p. But in general, sources of codimension p in gravity and electrostatics 
are well known to give singular solutions ~ r^^'°, so it should be no surprise to find a logarithmic singularity in the 
metric in the codimension-2 case. To study this in a controlled manner, we will regularize the singularity by giving 
the brane a finite thickness tq, and eventually take the limit rg — > 0. 

In terms of physics, our goal is to answer two specific questions: 

• Can the deficit angle dynamically relax to a value which cancels the effect of the brane tension? 

• Do these models lead to FRW cosmology when the matter on the brane has an arbitrary equation of state? 

The plan of the paper is as follows: in Sec. II, we review the argument behind the claim that one can only have pure 
tension on codimension-two delta function branes in Einstein gravity. In Sec. Ill we introduce the thick-brane setup 
we will be using to circumvent this difRculty. In Sec. IV, we review the known static background solution which we 
will be perturbing to obtain cosmological solutions. We present the linearized perturbative equations of motion and 
the Friedmann equations which follow from their solutions in Section V. Section VI rephrases these results in terms 
of effective four-dimensional quantities that are relevant to observers living on the brane, and derives the Friedmann 
equations. The reader who is interested in the main results of the paper may wish to skip directly to this point. We 
discuss our results in Sec. VII and give conclusions in Sec. VIII. Full solutions to the perturbed equations of motion 
are quite lengthy, and are provided to the interested reader in the appendices. 

It should be emphasized that supersymmetric versions of the model we are considering [l6l liM l4a | may be more 
promising, and represent a worthwhile subject for further research. 



II. THE TROUBLE WITH DELTA FUNCTION BRANES 

It has already been pointed out in 15] and confirmed by |2l| that in Einstein gravity, codimension-two delta function 
branes appear to be incompatible with the matter on the brane having a general equation of state. We repeat here 
the crux of the argument. 

We start with the most general six dimensional metric with the two extra dimensions having axial symmetry 



ds' 



N^{f, t)dP + A^{f, t)dx^ + B^{f, t)dr^ + C^{f, t)de^ + 2E{f, t)dfdi 



(1) 



where t and x correspond to the usual four large dimensions, while r and 9 span the two dimensional internal space. 
This metric can be put into the following form through a suitable transformation of the coordinates f[r^ t) and i(r, t): 



n^{r, t)dt^ + a^{r, t)dx^ + f{r, t){dr^ + r^d9^). 



(2) 



The equations relevant to our discussion are the (tt) and (xx) components of the Einstein equations, (here and in the 
rest of the text, primes stand for derivatives with respect to r while dots represent derivatives with respect to t), 
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where V^ is the two-dimensional Laplacian for the coordinates r and 6. Supposing there is a codimension-two delta 
function brane at the origin of the internal space, T^ and T^ contain the singular terms —p{t)5{r)/{2'iip{r,t)r) and 
p(t)(5(r)/(27r/^(r, t)r) respectively. Since in two dimensions the delta function can be written as 

J(r) = 27rrV2ln(r), (5) 

the singular parts of the Einstein tensor must come from the terms V^ ln(/), V^ ln(a) or V^ ln(n). It is not immediately 
obvious how to interpret having a(0,t) ~ r-'""^*^ or n(0,t) ~ r-''"^*) if we are to put the standard model on a 3-brane 
at the origin, since the metric is either vanishing or singular (depending on the sign of the exponent) there. If we 
consider this to be a problem, the Einstein tensor only allows for 

and we are led to conclude that codimension-two delta function branes can only have a stress energy tensor of the 
form Tj* = T^ (with T^ = Tg = 0). Furthermore, the deficit angle is related to an integral over the singular part of the 
Einstein tensor, and is therefore rigidly related to the brane tension through the Einstein equations. This prevents us 
from being able to meaningfully ask the question "what happens if there is a sudden phase transition on the brane, 
so that the tension and deficit angle become detuned?" unless we somehow regularize the behaviour of the metric at 
the brane. 

One possible regularization is to give the brane a finite thickness. By smoothing the singularity in the source, we 
remove the necessity of making the metric components vanish or become singular at the position of the brane. This 
is the approach we will follow in this paper, since we feel that it is less contrived than certain modifications to general 
relativity which could alternatively be invoked, and it allows for a better physical understanding of the apparent 
obstruction to zero-thickness branes with a general equation of state. 

Among the other possibilities that invoke modified gravity, which we will not pursue here, is the addition of an 
intrinsic curvature term to the brane action, as in Dvali-Gabadadze-Porrati braneworlds (see, e.g. pTl \4^). It has 
been argued that this scenario should arise naturally since the matter on the brane can induce a 4D kinetic term for 
gravitons through loop corrections 4^. Such terms can allow one to obtain additional delta function contributions 
in the equations of motion to match with a brane where p ^ —p. In this scenario, the deficit angle acts on the brane 
like a negative cosmological constant which, in contrast with what happens in the absence of the intrinsic curvature 
term, is no longer forced by the equations of motion to be equal to the brane tension. 

A second possibility, as was shown in |22|, is to add Gauss-Bonnet terms to the bulk action, which allows one to 
both put arbitrary types of matter and recover Einstein gravity on the brane. The result, as with the previous option, 
is that the exact cancellation between the brane tension and deficit angle is no longer imposed by the equations of 
motion. 

Recently a third proposal [40|,|4l| was made, using the intersection of codimension-one branes rather than a simple 
"point-like" codimension-two brane as the place where matter is localized. However since both references include 
extra terms in the gravitational action, it is not clear whether it is the different geometry or the generalized action 
which allows one to generalize the brane equation of state. 

Our view is that the above modifications are rather artificial because there is in general no reason to expect the 
metric to be well-behaved in the vicinity of a codimension-two or higher source. It appears to be an accident that the 
singularity is merely a conical one in the case of codimension-two source with pure tension. We will show below that 
the solution for a more general equation of state has more singular behaviour near the brane. 

III. SETUP 

As we have explained in the previous section, to allow for a general equation of state on a codimension-two brane 
with nonsingular solutions, it is necessary in Einstein gravity for the brane to have finite thickness. To study such a 
setup, we will model the branes using step function cores. This regularization scheme is similar to one that has been 
studied in the context of cosmic strings, (see, e.g. ^49,, 50j), with the difference that we will allow for the possibility of 
time dependent thickness. (This time dependence cannot be ruled out a priori for thick branes whose matter content 
itself is dynamical, as will be the case here). Furthermore, in order to find analytical solutions, we will treat the 
matter with p ^ —p as a perturbation to a static pure-tension brane background. We will therefore write the metric 
ansatz as 



with 



ds^ = -n{r, tfdt^ + a{r, tfdx^ + b{r, tfdr"^ + c{r, tfd0^ + 2E{r, t)drdt (7) 

n{r,t) = e^«('')+^i(''^*^ a(r, t) = ao(Oe'^"^'^'+^''''^*^ 6(r, t) = 1 + Bi(r, i); 

c{r,t) = e^o^'-'+'^i^'^'*); E{r,t) ^ Ei{r,t). (8) 



The bulk stress-energy content consists of a 6D cosmological constant Ag and a two- form field Fab, so that the bulk 
action can be written as 57] 

^ - As - ^F'^'^Fab) . (9) 



Although there exist solutions for both positive and negative values of Ag |la|, we will take it to be positive only 
for the rest of the paper. We assume that the only nonvanishing component of the vector potential is Ag{r,t) — 

Ag '(r) -I- Ag {r,t). For the sake of generality, we include a possible perturbation of the 6D cosmological constant, 
Ag ^ Ag -I- <5Ag . The background flux from the gauge field is necessary for getting a simple spherical compactification 
of the two extra dimensions. The full stress-energy tensor is taken to be of the form 

T,^ir,t) = ttir,t) + 9iro{t)~r)SSir,t) + 9ir-r,{t))S4ir,t). (10) 

Here, ro{t) and r*(t) are the radial positions of the edges of the thick 3-branes centered around the poles of the 
compact internal space, i^ refers to the bulk content, and S^ is the core stress energy, given by 

Si = -(T-p; S'^ = -a+p; S'^^O + p^; 

S'g = 0+p'g; Sl^O+pl; S', = + pU 

S4 = -a - p*; S^l ^ -a + p^ S'*': = + p4; 

S4 = 0+p4; S4 = 0+p4: 5,*-0 + p4. (11) 

We treat the time dependence of the thickness as a perturbation, so that ro(t) = ro -|- Aro(i), r*(t) — r^, + Ar*(t) and 

0(ro(t)-r) = 0{ro - r) + Siro ~ r)Aroit) + 0{Arl) (12) 

e{r-r^{t)) = 0{r~r^)~6{r^r^)Ar^(t) + O{Arl). (13) 

so that effectively, we can write the stress-energy tensor as 

t1: + s1; + s4 (14) 



with, e.g., 



= -a9{ro-r) + [-pe{ro-r)-a6{ro-r)Aro{t)] (15) 

= ~ae{ro~r) + [pe{ro-r)-a6iro-r)Aro{t)] (16) 



and similarly for all other terms. 

Here a represents the tension of the regularized brane, and p,p represent contributions from ordinary matter on 
the standard-model (SM) brane, while starred quantities refer to matter on a hidden branc which is antipodal to 
the SM brane on the two-sphere bulk. All variables except a are functions of both r and t. The appearance of 
off-diagonal terms in the perturbative stress-energy tensor might be surprising at first, but as we will explain, they 
must be included, since they cannot be made to vanish in a coordinate invariant manner. 

The subscripts on the metric and gauge field perturbations indicate their order in a perturbative series in powers 
of p.|58| We will furthermore assume that time derivatives of the perturbations are of 0{p^^^), which is implied by 
the usual law for conservation of energy p ^ (a/a)p ^ p^'^ . The framework is thus the same one that was used in the 
context of Randall-Sundrum cosmology in [52.l53| . 

IV. THE BACKGROUND 

In this section we will generalize the static "football shaped" solution jla . ll4l | to the case of branes with nonvanishing 
thickness. (This will serve as the background around which we will later perturb by adding cosmological matter on the 
branes.) We will show that in the limit where the brane thickness vanishes, we recover the solution for delta function 
branes. This might at first seem like an obvious outcome, but in fact it does not necessarily happen, in the context 
of codimension-two objects in Einstein gravity |54| . In the present model, the extra dimensions are compactified 
on a two-sphere, with pure-tension branes positioned at the poles. The only effect of the branes on the bulk is to 
induce a deficit angle so that the internal space effectively looks like a sphere with a wedge taken out, hence the name 
"football-shaped" (the European reader may prefer to think in terms of a rugby ball). We will refer to the interior 
region of the thick brane as the "core," and the exterior as the bulk. 



A. Core solutions 

We first consider static solutions of the unperturbed Einstein and gauge field equations 

SirGeia + Ag) + inGee-^^^A'-^^'^ + SA'^ + C'^ + QA'i^ + C^ + 3A^C^ = (17) 

< - Al + < _ 3^^2 ^ c-^^^, _ ^, ) ^ 2^^7V, ^ (^8) 

STrGeAe - A^G^e-^^" Af^'^ + 3A^(A;, + A^o + Cq) + Cq^o = (19) 

3A[; + 7V^' + 3A^2 _^ ^/2 _ ^/^/ _ 3^/^^/ ^ (20) 

a(3^;, + A^^) = (21) 

Af + Af[ZJ^, + N',-C',) = Q (22) 

in the interior (core) of the thick brane centered at r = 0. In the absence of any matter perturbation, conservation 
of energy implies eq. H21|) . Once this condition is imposed, combining the [rr) and [69) Einstein equations (|18|l and 
(|20|l gives ^Q = A'q = 0, and the system of equations is solved by 

gCoW ^ isin(fcr); 
k 



with 



Mr) 


= 0; No{r)=Q: 


Cir) 


= A(cos(fcr)-l) 




k" =%ttGq{<7 + P'^): 



(23) 



P^ = 2A6. (24) 



In addition to this static solution, there exist solutions where the brane worldvolume is deSitter or anti-deSitter 
space if we relax the tuning between the gauge field and cosmological constant [l3- In the context of a related 
supersymmetric model, it has been shown |l(l l45l | that the flat brane solution is actually singled out by the dilaton 
equation of motion. Nothing forces us to choose the flat brane solution, but we do so for simplicity. It should also be 
noted that our background assumes both branes have the same tension. While there do exist solutions with different 
tensions on the two branes, these include warping and cannot be put in closed form, so the choice to expand around 
the simple unwarped background is also made for simplicity. Since our framework will allow us to perturb these two 
unnecessary but simplifying tunings (flat branes and equal tensions), our final results will not be affected by this 
choice. 

We assume that the system is symmetric around the equator of the extra dimensions. Thus there is a brane also at 
the south pole, whose coordinate position is taken to be r = t: /k — <j). The solution in the core for the second brane is 

Ao(r) = 0; NQ{r)=Q- e^«W = ^ sin(fc(r + 0)); 

k 

4°^W = ^ (cos(Mr + 0)) + 1) (25) 

where the difference in sign relative to the solution in the upper hemisphere is necessary for the gauge field to vanish 
at both poles. The phase (p is an integration constant to be determined by matching the core solutions smoothly to 
the bulk solutions, which we now consider. 

B. Bulk solutions 

In the bulk, the static solutions look like those for core,|52 except for the difference that a — Q. The bulk solutions 
are therefore 

gCoM ^ J_ gin(^(^ + 0)) (26) 

ki 

aT {r) = /3^ cos(fc(r + 4>)) - Pi (27) 

Demanding that the metric and gauge field as well as their first derivatives be continuous across the junctions located 
at r = To and r = 7r/fc — tq — (j) = r^ (assuming that both branes have the same thickness), and using the fact that 

P^ =p^ = 2Kq- P = SnGeP^ (28) 



(c.f. eq. (|24|l ). we find the following matching conditions from the r = tq junction: 

j-_2 _ sin(/i;ro)^ cos(fcro)^ 



1 _ fc2 

- - A: 

tan(fc(ro + (/))) = — tan(fcro) 



fc2 



h 



A:2 



1 — cos(fcro) 1 — 



U2\l 



fc2 



(29) 
(30) 
(31) 



In the southern hemisphere, the bulk solution has the same form as for the northern one, except that a different 
constant of integration (31 must be used in place of /3i in the solution for the gauge field. At the r — r^, junction we 
obtain 



~K' = 



sin(fcro)^ cos(fcro)^ 



/c2 



B 



' k - - - -\ k 

tan I — TT — Zero — fc0 + fc(/) J = — — tan(fcro) 



/3* 



k 

/fc2 



1 — cos(fcro) 1 



(32) 
(33) 

(34) 



The matching of the bulk gauge field solutions to both cores requires that f3l = — /5i, which makes the gauge field 
discontinuous at the equator. If the gauge field couples to matter, such a discontinuity is physically sensible only if 
the solutions in the two hemispheres are related to each other by a single- valued gauge transformation at the equator, 
which leads to the quantization condition (see p!a.l55|') 



Pi 



2ff 



where n is an integer and g is the U{1) coupling constant. 
Finally, we can use the relations (|30|l and H33|) to show that 



- TT mn 



(35) 



(36) 



where m is in principle an arbitrary integer, but the choice m = — 1 ensures that both and 4> vanish when a vanishes. 
Let us now show that the usual relation between the deficit angle and the brane tension emerges in the tq — *■ limit 
of the above solution. Since kro remains finite, the brane tension a scales like 1/rg. The effective four-dimensional 
tension is found by integrating over the volume of the extra dimensions (see Sec. VI): 



.(4) 



27rcr 



(1 - cos(fcro)) 



so that as ro ^ 0, 



AGea^^^ = (1 - cos(fcro)). 



(37) 



(38) 



We must now consider how the deficit angle should be defined in the case we are considering, where the brane is thick. 
From the bulk point of view, the radial distance from the brane at r = R is R — ro. The circumference of a circle of 
radius R is 27rc(i?, i), while the circumference of the brane is 2'Kc{ro,t). If there is no matter on the brane, so that 
the internal space is perfectly spherical, we would expect that as tq ^ and i? — > 0, 



27T[c{R,t)-c{ro,t)]^2TT{R-ro). 

On the other hand, if there is matter on the brane, it will modify the previous relation to read 

A 
2^ 



2^(c(i?, t) ~ c(ro, 0) = 2^(E - ro) 1 



Thus we can define the deficit angle as 



A EEE 27r lim 



lim 1 



c{R,t) -c(ro,i) 



R-ro 



(39) 



(40) 



(41) 



To lowest order, this simply leads to 



27r ( 1 — lim lim 



fcf^sin(fc(i? + 0)) 
R 



27r 1 - lim 
2tt {1 — cos{kro)) 



cos(fcro) sin(fci?) 



kR 



The expected relation between the deficit angle and brane tension follows: 



(42) 

(43) 
(44) 

(45) 



We have therefore shown that this thick brane model consistently yields known results for a delta function branc in 
the limit of vanishing thickness. In other models of localized stress energy, it could happen that the (rr) and {96) 
components of the stress-energy tensor do not vanish in this limit, in contrast to the stress tensor for a brane. In 
such a case, the relation (|45|l would no longer be guaranteed |5J| . Below we will demonstrate another means by which 
the relation (|45|l can be relaxed: adding brane matter with a different equation of state generalizes the form of H45|l . 
This is derived for thick branes, and it will be shown that the thin-branc limit is singular for these solutions. This is 
probably why departures from l|45|l have not been pointed out before, to our knowledge. 



PERTURBATIONS 



It would be very interesting if an inflating solution, which started out with some relation between the brane tension 
and deficit angle differing from that required for a static solution, would dynamically relax to the static solution. The 
demonstration of such behaviour would be a convincing step toward solving the cosmological constant problem. It 
would also suggest that cosmology should be modified for ordinary matter on the brane — given that there is a drastic 
modification when the source on the brane corresponds to vacuum energy. To explore these issues, we want to add 
possibly time-dependent perturbations p and p to the stress-energy of the brane and see how the geometry responds. 
A physical example where this would be relevant is a phase transition on the brane, where part of the tension gets 
converted instantaneously into radiation. 

Before considering the perturbed equations of motion, it is useful to find a set of coordinate invariant quantities 
associated with the metric perturbations in ^ . Considering a small reparametrization of r and t such that 



r -^ r + f{r,ty, t^t + h{r,t), 

where / and h are assumed to be of 0{p), the perturbations transform as|6(]| 

iVi -> A^i; Ai^Ai; Bi^Bi+f; Ci -^ Ci + C'J; 



(46) 
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p\ - ah'. 







s\ -f a5{rQ - r)f; 



Pt 



Pl + 'yf; 



(47) 



One can see why it is necessary to include off-diagonal terms in the perturbations to the core stress-energy tensor S^, 
: they cannot be made to vanish in every frame. Moreover to have a consistent perturbative series in powers of p, 
_Ei,p[,p* and h are C'(p'^/^), while all other terms are of order p. 

From the form of the gauge transformations, we can deduce the following 11 gauge-invariant variables: |6l| 
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= N'l-A'i; 




W = 3A'i + N'l 




r-il (-Ml 




Y = A^y 


~Af\Bi^Ci)- 
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= p + S{ro - 


r) 


—Ci+aAroit) 
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p = p- (5(ro - r) 


—Ci+aAroit) 
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r Ci 

- Pt f^^,; 




pI=pI+ pI - a 


El. 


P5=Pr; P6 


-Pl-fr- 



(48) 



In terms of these variables, the equations of motion at first nontrivial order in the perturbation are 62 



Z' + C'oZ = 
W'-CqW = 



p', - Cl^Pe + aW 
U' -Y - (3e^"X 

ao 4 



SirGpe 




+ 87tG{p + p) 



Oo 
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+ 87rG(p5 - SAe) 



-C'^W - 8TTGPe-^°Y 
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-27rG(/5-3p + 3^6 + 4^6) 



SirGiPl - (3e-'-'>U) 



Pt 



C'oPl + oX 



(49) 

(50) 
(51) 

(52) 



(53) 

(54) 
(55) 
(56) 

(57) 
(58) 



These are not all independent however. It can be shown that cq. H51|l and eq. H55(l can be obtained as combinations of 
the other equations in the system, and can therefore be ignored. We have found the general solution to this reduced 
system of equations in appendix A. 

The general solution to (|49I58|I is complicated, but we are principally interested in just one aspect, the Friedmann 
equations which describe the expansion of the universe for an observer on one of the branes. In the present construction, 
the Friedmann equations arise through imposing the appropriate boundary conditions at the junctions of the regions 
interior and exterior to the branes. In terms of our gauge-invariant variables, continuity of the functions but not of the 
derivatives is required at the interfaces, since the equations of motion do not contain any spatial second derivatives. 
Moreover we require that all functions be nonsingular at the poles. Finally, in what follows, we assume that p,p, /5,, 
p* are functions of time only and that pe = e^^°'^''^Ve,{t) and p*6 = e^^°^''^V^&{t).^ 

With this ansatz for the brane energy density, applying boundary conditions allows us to derive the Friedmann 
equations. However, there is one last point to address before writing them down: the 6D quantities p and p we have 
used so far are not the ones an observer on the brane would identify as the energy density and pressure. It is therefore 
necessary to define what the effective energy density and pressure are for a 4D observer if we want to compare the 
Friedmann equations in this model to the standard ones. 



VI. EFFECTIVE FOUR DIMENSIONAL QUANTITIES 

To define the effective four dimensional quantities, we integrate the 6D quantities over the thickness of the brane, 

/■ra 

S^^'^l = 27r / dr b{r, t) c{r, t) S^^^l (59) 



which perturbatively leads to 



a(4)+p(4) = 27r / dre^°{l + Bi+Gi){-sl) 
Jo 

pro 

_cr(4)+p(4) = 27r / dr e^^^l + Bi + Ci){sl) 
Jq 

/•ro 

pW ^ 2tt dre^°p5 

Jo 

pro 

p^^> = 27r / dre^"p6. 
Jq 



(60) 
(61) 
(62) 
(63) 



and similarly for the other brane.|6J| In order to compute p^^^ and p^^\ we choose a gauge and find Ci or Bi 
explicitly. For this purpose it is convenient to work in coordinates where Bi — Ei =0. From the form of the gauge 
transformations (|47|l . it can be shown that the 4D quantities are gauge invariant. For example, 5{(7'^^^ + p'-"'-') = 
27r JT" dr -^{e'~'°af) is a surface term which vanishes because / = at r = and a — at r = r^ . 
Similarly, the 4D Newton constant is related to the 6D one by dimensional reduction, 

Ge = GixV (64) 

/'27r p7r/k — (p 

= Gi dO e"^°('') (65) 

Jo Jo 

Air - - 

= G4 X 7^T^(fc' + (fc' - P) cos(fcro)) (66) 

where we neglect corrections of 0{p). 

After gauge-fixing, the Fricdmann equations and conservation of energy can be rewritten in terms of the effective 
4D quantities (see the appendix for the full solutions) 

^)' = ^(p^'^+4'^+Acff) (67) 

^ = f^y_4.G4(p(^)+p(^)+pl^)+pl^)) (68) 

p(4) ^ _3^(p(4)+p(4)A (69) 

pl^^ = -3^(pl^)+pl^)). (70) 

These constitute the main result of this paper, and show that we indeed recover standard cosmology on the brane, 
except for the constant of integration Aoft which represents bulk contributions to the 4D cosmological constant, and 
the contribution of "dark matter" which is on the second brane. Of course we also expect corrections of order p^ 
and higher which would appear at higher order in the perturbation expansion |53j| . but these are not relevant for the 
immediate question concerning self-tuning. 

VII. DISCUSSION 

To elucidate our results H67I70|I . let us first explain the role of the constant Aeff. It might at first be surprising 
that our solutions contains a term that does not vanish when matter perturbations are absent. However, this simply 
reflects the fact that our choice to perturb around a static background solution was arbitrary, and there exist evolving 
solutions even in the absence of matter perturbations. This is because we are free to choose any set of values for the 
parameters of the background solution a'''^\ (3 and Ag.JGa] We arbitrarily tuned 0^ = 2A6 to ensure that the brane 
does not expand. Solutions with nonzero Aeff simply correspond to using a different set of values for the background 
parameters, for which the tuning between the field strength and bulk cosmological constant has been perturbed. The 
value of the constant K^s is thus determined by the background around which matter perturbations are being added, 
and cannot be used to subsequently tune the resulting solutions. 

We will now consider some limiting cases which provide a consistency check on our solutions, and which illustrate 
the differences between the thick-brane model and its thin-brane limit. 

A. Thin brane limit 

First we would like to consider how our results behave as we take the thin brane (i.e., tq -^ 0) limit. Earlier we 
argued that for a delta function brane with an equation of state differing from pure tension, the 4D components of 
the metric have singular behaviour at the origin. Let us explain how this is manifested in our results. 

The thin brane limit corresponds to taking k — > cxd, since the constant kro remains finite. This follows from the 
requirement that the effective four dimensional quantities cr^ \ p^'^' and p^'^' should be held fixed as we vary the brane 
thickness. We can then see from our results why the ro ^ limit produces singular results at the branes unless the 
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matter on the brane consists of pure tension. From the metric perturbations Ai and A^i, we find that|66 



N[iro,t) = 3fc^A^--^+0(l/fc) (72) 

sm(fcro) 

where we have concentrated only on the brane located at the upper pole, but would of course find similar results 
for the other brane. These diverge in the thin brane limit (since k '^ l/j'o ^s rg — > 0) unless the perturbations 
satisfy p^'^' = —p^^'. Thus it is only possible to have pure tension on a codimension-two delta function brane, if we 
(unreasonably) require that gtt and gxx be regular at the brane. On the other hand if we allow an arbitrary equation 
of state, then A[ and N[ diverge like l/r near the brane, indicating that Ai and iVi go like In r. 

We stress again that regularity at the brane is probably too much to ask of a solution when a delta function source 
is present; rather one should expect singularities in this case. Therefore, we do not interpret our results as indicating 
any fundamental obstruction to having general types of matter on a codimension-two delta function brane. Rather, 
eqs. (|71I72|) imply the mildly singular behaviour 

a{r, t) - r-^'^(p+P'> , n{r, t) - r^-^^ef^+P) (73) 

for the 4D metric components in the vicinity of the brane. 

B. Deficit angle 



We turn our attention now to the deficit angle which, as was previously mentioned, is defined by 

c{R,t) - c{rQ,t)' 



2iT lim 



lim 1 



(74) 



70^0 R-ra 

around the brane located at r = 0. Plugging our solutions c{r, t) « e'-^°('')(l + Ci(r, t)) into this, we obtain 

A = 2^G6(4fT(4)+p(4)-3p(4)) (75) 

The same procedure performed around the other brane leads to 

A, = 27rG6(4a(4)+p(4)_3p(4)) (76) 

We see that the deficit angle indeed responds to changes in the stress-energy on the brane; moreover in the familiar case 
where p^^^ — ~p^^^ (pure tension), the relation reduces to its expected form A = SnGeaj-gl where alJ^ = a^"^^ + p^'^K 
However, this is not enough to insure a static solution, as can be seen from eq. (|67|l . At first this looks quite mysterious: 
in the pure tension case, why is the effect of the tension not canceled by the conical singularity, as was the case for 
the unperturbed solution? The answer, given in the next subsection, has to do with the effect of the perturbation on 
the bulk fields, in particular the gauge field. 

C. General equations of state 

With our results we can now answer the physical question posed at the beginning of section 5: suppose we start from 
a static solution which undergoes a phase transition that converts a part of the brane's tension into radiation — how 
does the geometry respond to this sudden change in the equation of state? The perturbation to the stress tensor can 
be thought of as the sum of two pieces, 

p = Sa + piad 

P = -5a + -prad (77) 

where 5a is negative, and initially p ~ Q. According to the Friedmann equations, the universe will start to collapse 
toward a big crunch, as one would expect for an anti-deSitter solution; the universe behaves just as though it had a 
negative cosmological constant. The conclusion actually does not require the transition to be sudden, nor does it have 
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to produce radiation. The same outcome occurs even if the equation of state of the perturbation changes arbitrarily 
slowly from w = —1 to any larger value. 

For example we could consider such a transition in a more smooth way by choosing a positive initial value of dp, 
tuning the solution (using the integration constant Aoff ) to be initially static, and imposing a time-dependent equation 
of state for the perturbation, w — (— 1 + tanhet)/2. This represents the gradual conversion of the perturbation from 
vacuum energy in the past to pressureless matter in the future. Again, the Friedmann equations predict a collapsing 
universe due to the effectively negative cosmological constant at late times. 

Similarly, if we start from an initially static solution with a negative value of 6p and the opposite time dependence 
for the equation of state, w = (—1 — tanhet)/2, we obtain an inflating solution at late times due to the positive change 
in the brane tension. In all cases, it is clear that no self-tuning mechanism cancels the effect of the change in tension. 

Now we address the seeming paradox: how can our results be self-consistent? We find that the Hubble rate is 
sensitive to changes in the energy density on the brane, yet the deficit angle tracks this energy density in the way 
which we would hope for if there was self-tuning of the 4D cosmological constant to zero. The key point is that the 
tuning between the deficit angle and brane tension is not sufficient to lead to a static solution; one also needs the 
field strength and bulk cosmological constant to be tuned. In general, Sp induces a perturbation to the field strength 
F^ = Fq + 2Fq5F , where in terms of our variables, we find that 

<5F = \/2e"^«r(r,i) (78) 

In the static solution, the field strength had to be tuned relative to Ag. Apparently it is this tuning which is spoiled 
in the presence of matter on the brane, rather than anything related to the deficit angle. 

VIII. SUMMARY AND OUTLOOK 

In this paper we constructed a thick codimension-two braneworld model in order to determine whether such scenarios 
have a self-tuning mechanism which could help solve the cosmological constant problem, and also to investigate the 
question of whether they allow us to recover standard cosmology. 

The good news is that the second question is answered in the affirmative. Codimension-two braneworlds are quite 
consistent with the general expectations from decoupling and dimensional reduction: the 4D effective theory looks like 
general relativity at lowest order in the density of matter on the branes. Although corrections to GR will undoubtedly 
appear at order p^, we have not tried to compute these in the present work. 

On the other hand, we find an obstacle to self-tuning of the effective 4D cosmological constant in this model. Sarting 
from a static solution, any phase transition which changes the tension of the brane will either lead to a contracting or 
an expanding universe. In the second case, once the matter has redshifted away, we do not recover the original static 
solution, but rather a deSitter one, where the difference relative to the original static solution comes from the fact that 
the phase transition has spoiled the original tuning between the bulk gauge field strength and the bulk cosmological 
constant. 

We were motivated to study a thick brane model because of the apparent impossibility of putting matter with a 
general equation of state on a codimension-two brane in Einstein gravity. We find that there is actually no prohibition; 
instead the metric becomes singular near the brane when it contains general kinds of matter. In retrospect this is not 
surprising, and the fact that pure-tension branes induce only a conical singularity appears to be fortuitous. Accepting 
the singular nature of the metric, we derived an expression for the deficit angle with arbitrary matter on the brane, 

A = 27rG6 f4a(4) + p'^) _ 3^(4) A (79) 



which reduces to the expected one for a pure-tension brane. 

Although our results show that codimension-two braneworlds in Einstein gravity don't solve the cosmological 
constant problem, the outlook may be better in models where there is supersymmetry in the bulk [I^ |43, |4a| . The 
encouraging point in our results is that the deficit angle does respond to arbitrary changes in the brane stress energy 
in such a way as to cancel its contribution (through the singular part of the gravitational action). What is needed is a 
symmetry to maintain a vanishing bulk contribution to the effective 4D cosmological constant, and it has been argued 
that supersymmetry can do precisely this. In SUSY models, the dilaton equation of motion insures that the tuning 
between the gauge field and bulk cosmological constant (which appears in the guise of a scalar field potential) required 
for a static solution always holds as long as the dilaton is stabilized. Since we have shown that it is the detuning 
between the field strength and bulk cosmological constant that leads to non-static solutions in Einstein gravity, it 
seems reasonable to expect that in models with supersymmetry, there will indeed be a self-tuning mechanism to cancel 
the effective cosmological constant. Checking this explicitly is the subject of ongoing work. 
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APPENDIX A: PERTURBATIVE SOLUTIONS 



We present here general solutions to the perturbed equations of motion (|49I58I) . The superscripts /,//,/// and 
IV are meant to identify quantites associated with the different regions in the model, namely the upper core (r ~ 
to r = ro), the upper bulk (r = rp to r = 7r/2/fc — </)), lower bulk (r = Tr/2/k — (j) to r — 7r/fc ~ cj) ~ tq) and lower core 
(r = ir/k — — ro to r = n/k — (p). We have defined f = r + (j) and f = r -\- (j). We have not written the solutions for 
the gauge invariant variable U{r,t), since it can be obtained trivially from the algebraic equation l(57|l once the other 
variables have been solved for. 



Z' - 



Y' 



X' 



Pt 



Tiit) _^cos{kr) 



sin(fcr) ksa\{kr) 



47rG6(p + p) 






W' = J^2^(i)sin(fcr) 



STrGg sin(fcr) cos(fcr) 



tIi^-\ , '^cos(fcr)^ (P + 87rG6cr)(cos(/i:r)2 -sin(fcr)2) 
P5 - -/"a v-) "f z -^2 W -tta ^6 



4fc4 
(fc2 + STrGgcr + 2{P - SirGecr) cos{kr)'^) sm{kr) 

4fc5/3 
(fc^ — SirGea) sin(fcr) cos(fcr) 



sin(fcr) 
SnGePk 

Hit) ^ 
cos(fcr)2 

irGe 
I6k^ cos{kr 



SnGepk^ 



Hit) 



flo 



3— + 3 — + SirGeHit) - STrGgMg 
ao \ao 



(5fc^ — IGnGea) cos{kr] 



6fc3 



Hit) 



[I7fc2 + mnGea + 32 cos{krfik^ - inGea) 



8cos(fcr)4(llfc2 - leTrGgcr)] Ve 



{2cos{krf -1) 



8fc2cos(fcr)2 



9— +9 — -4^G6(p-3p + 8M6-4.fJ(i)) 
ao V«o 



cos(fcr) (iQ 
~3 . . —ip + p) 
ksva[kr) ao 

cr7rG6-2cos(fcr)2(fc2 



d_ 
di 



9cr(l + 2cos(fcr)2; 

8k^ sin(A;r) cos(fcr) \ ao 



airGe) airGeil + 2cos{kr 



2k^ sm{kr) cos(fcr 
crvrGg 



k^ sin(fcr) cos(fcr 




2Hit) 



ASk"^ sin(fcr) cos{kr 

8 cos(fcr)4 (llfc2 - leanGe)} Vq 

(5fc2-167rG6cr)crcos(fcr)2 



{-3(17fc2 + IGcnrGe) + 96cosikrf{k^ - 4cr7rG6) 

Hit) 



12fc4 sin(fcr) 



-Hit) 



ksm(kr) cos(fcr) 



Hit) 



sin(fcr) 



(Al) 

(A2) 
(A3) 



(A4) 



(A5) 



(A6) 



7ll 



Ti\t) cos(fcf) 



sin(A;r) fcsin(/cr) 



ao 
ao 



(A7) 






II fcsin(fcr) cos(fcr) jj , , sin(A:r) 



SnGePki 



TT 5 cos( kr) IT , , 



X" = 



+ 



6fc 
cos(fcf)2 



87rG6/3fci 

(2cos(fcf)2- 1) 
8P cos(fcf)2 oo 






9^ + 9^"" 
ao 



- 327rG6(5A6 
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(A8) 
(A9) 



(AlO) 



Ziii 

ylll 



H'\t) ^ cos(fcr) 



sin(/cr) A:sin(fcr) 

J^2"^(t)sin(fcf) 

k sin(fcf) cos(fcf ) jjj 



Oo 

ao 



X 



III 



6fc 
cos(fcf)2 



STrGe/^/ci 



•^2^^nt) + 



sin(A:r) 



5cos(A:r)^„^ 



STrGg/J/ci 

(2cos(fcf)2- 1) 
8fc^ cos(fcf)^ 



ao 



ao 



3^ + 3f^) -STTGeMe 



- 327rG6<5A6 



(All) 
(A12) 
(A13) 



(AM) 



^iv ^ 



Ti^{t) cos(fcf) 



-2 



sin(/cf) A:sin(fcf) 



4.G6(p*+p.) + ^-P 
ao V ao 



P*5 
ylV 



tIVu\ ■ fi '\ SttGo sin(A:f) cos(fcf) 
J^2 (t) sm(fcr) p V*6 

tIV(.\ , <7Cos{kf) wy/,N (fc^ + 87rG6cr)(cos(fcf)2 - sin(fcf)2) 
■^3 W + — y-2 (t) ^4 '*e 

(fc2 + STrGgCT + 2{k^ - SirGea) cos(fcf)2) sin(/cf) 



4fc5/3 
(fc^ — STrGecr) sin(A:f) cos{kf) jy 



P*e 



sin(fcf) 
SirGePk 



SirGelSk^ 
3^ + 3^"° 



•^r(i) 



ao 



ao 



SirGeTsitJ-SnGeSAe 



X 



IV 



Ti^'it) {5P-16TTGea)cos{kf)_jy 



6fc3 



cos(fcf)2 

16fc^ cos(fcf)2 
8cos{krf{nk'^ - IGTrGgcr)] "P, 



•^2 W 



[I7k^ + 16nGe<7 + 32 cos{krf{k^ - inGecr) 



~ r 

p*t 



(2cos(fcf)2- 1) 
8fc2cos(fcf)2 

cos(fcf) do 



,00 

ao 



9— + 9 ( — ) - 47rG6(/5, - 3p, + 8(5A6 - 4jc-j^(f)) 



= -3 "77,' U —if* +P*) + 77 
/esin(/er) ao ai 



9(T(l + 2cos(fcf)2) I do 



ao 



8^;"^ sin(fcf) cos(fcf) \ ao \ao 
a7rGe-2cos{krf{k^ -cnrGe) crTrGgll + 2cos(fcf)2) /3 



2fc3 sin(A:f) cos(fcf) 
cTTrGg |- 0/1 f7j 2 



■p* 



p, + 2^3^^(t) 



48fc^ sin(fcf) cos{kf) 

8 cos(fcf)''(llfc2 - IGcTTrGe)} "P, 



k^ sin(fcf) cos(fcf) \^ 2 
{-3(17fc2 + IGcTTrGe) + 96cos(A:f)2(/c2 _ AaTrGe) 



(A15) 

(A16) 
(A17) 



(A18) 



(A19) 
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(5fc — 167rG60-)cr cos(fcf 



f.\2 



•^rw 



-T^i"(i) 



sin(fcf) 



(A20) 



12k^ sm{kf) k sin(fcf) cos(fcf) 

These solutions solve the equations of motion H49I58|I . However, we still need to impose boundary conditions. 

Before doing so, we will specialize to a gauge where Ei = Bi = 0. We can then use the above results to find Ci(r, t) 

and Ag (r, t). (The explicit solutions for Ai{r, t) and Ni{r, t) are of no particular interest, so we won't bother writing 
them explicitly). 



ci = Hit) 



kcos{kr) 
sin{kr) 



Hit) 



(5fc^ — lecTTrGe) cos(fcr) ^ 



6fc3 



Hit) 



ttGr 



[sin(/cr)(17fc2 + IGcnrGe ~ Acosikrf{nP - IQanGe)) 



16fc^ sin(fcr) 

4(3fc2 + 16anGe)krcos{kr)] Vg 

2kr cot{kr) — 1 



8fc2 



9^ + 9 Jl _ 4^Geip - 3p + SMg - 4.F|(i)) 



Ci^ 



II 



H'it) + 



kcos{kr) 
sin(fcf) 



•^6 W + 



2krcot{kr) — 1 



8fc2 



5 cos(fcr) 
6fc 

2 



•^2 W 



9^ + 9 (^) -32^G6<5A6 



2kr cot{kf) 



kcos{kr) 
sm(kf) 



H'\t) 



5cos{kr 



H''it) 



8fc2 



ao 



- 327rG6M6 



Gr = ^i"(i) + 
ttGq 



•^2^^(i) 



+ 



fccos(fcf) _/v'. (5fc2 — IGcTTrGe) cos(fcf) 

sin(fcf) ^ ^^'^ 6fc3 

[sin(A:f)(17A:2 + IGcnrGe - 4 cos{kf)^ {Uk^ - 16a7rGe)) 
llfc^ - 16(jTTGe)kc/)) cos(/cf)] P^g 

IV 



16fcS sin(fcf) 

4((3fc2 + 16(T7rG6)A:r 

2fcrcot(fcf) — 1 



8fc2 



9— + 9 I — ) - 47rG6(p* - 3p* + SSAq - 'iTi'' (t)) 



(A21) 



(A22) 



(A23) 



(A24) 



.(1)^ ^1 , X /3sin(fcr) ^r, , Bcos(kr) ^r , , 

4^ = Hit) - ^^ ^ .^e'w + j^; ' Hit) 

{3k^{k^ ^ SanGe) + 'jp^nGej^k^ - legTrGg)) cos(fcr)2 
''' ASnGePk^ 

(4/327rG6fcrsin(fcr) - cos(fcr)(fc2 - G/S^TrGe)) 



Hit) 



k^(3 



Hit) 



1 



[(-8fc2(fc2 - 8cr7rG6) - 4/327rG6(llfc2 - IQanGe)) cos(fcr)^ 



48A:8/3 

(12fc2(fc2 + ScTTrGe) - i0^'KG(,{2%k^ + 80cr7rG6))) cos(fcr) 

127rG6/32(3fc2 + 16cr7rG6)fcrsin(fcr)] ^(0 

3(6/327rG6fcrsin(fcr) - cos{kr){P - 9(3^t:G6)) [ao Ao 
Sk^PnGe ao \ao 

l3TTGei2kr sm{kr) +3cos(fcr)) 

2fc4 
8f3^TTGekrsm{kr) - cos{kr){k^ - UP^nGg) 



ip-3p) 



k*(3 



Mfi 



i(i 



II 



= H'it) 



/3sin(fcr) 



H'it) + 



/3cos(fcr) 
kki 



H'it) 



(A25) 



15 



A 



(1) 



III 



i(i) 



IV 



-3P + (3fc^ + 20^^7rG'6) cos(fcf)^ 
3{6P'^TTGekrsm{kf) - cos(fcf)(P 



- 9f3^TrG6)) 



8fc3/ci/37rG6 
8/3^7rG6fcrsin(A:f) — cos(fcf)(fc^ 






UP^ttGg) 



Tr{t) 



(3 sin(fcr 



k^kiP 



SA, 



(3coa{kr) 



fcfci 



H''it) 



-3fc2 + (3fc2 + 20P^ttG6) cosikff 



■H"{t) 



ASirGePk^ki 

3(6/3^7rG'6fcrsin(fcf) - cos(fcf)(P - Q/J^TrGe)) 

8k3kif3TTGe 
8(3^TTG6krsm{kf) - cos{kf){k^ - l2f3^TrGe) 



Oo 



= ^rw- 



k^kifi 

/3sin(fcf) jy 



Mfi 



•^rw + 



7 KN I -^6 

(3fc2(fc2_8a7rG6)tan(fcf)2 



- 4/327rG6(5/c2 - IGaTrGg)) cos(fcf)2 



487rG6/3fc5 
(4/327rG6fcrsin(fcf) - cos(fcf)(fc2 _ G/S^TrGe)) jy 



^rw 



/(:4/3 



•^rw 



1 



[(-8fc2(fc2 - ScTTrGg) - 4/327rG6(llfc^ - lecrvrGe)) cos(fcf)3 



48A:8/3 

- (12A:2(fc2 + 8cr7rG6) - il3^T:G(i{29k^ + SOcrTrGe))) cos(fcf) 

- 127rG6/32fc(r(-3fc2 - lecrvrGe) + ^(-llfc^ + lecrTrGg)) sin(fcf)] V^eit) 

2' 



3{6P^TTG6krsm{kr) - cos{kf){k' - Q/g^TrGg)) 

8fc4/37rG6 
/37rG6(2fcrsin(fcf) +3cos(fcf)) 



Oo 
ao 



2fc4 
8/327rG6fcrsin(/cf) 



cos(fcf)(fc2-12/327rG6) 



A:4/3 



5K 



(A26) 



(A27) 



(A28) 



The conditions we will impose on the solutions are the following: they must all be regular at the poles, and must 
satisfy the jump conditions at the boundaries between the cores and the bulk. Furthermore, all solutions must match 
smoothly across the equator separating the upper and lower core, with the exception of the gauge field perturbations, 
which as explained in the main text, must be related by a single valued gauge transformation. 

The important point to watch out for comes from the equation of motion for the gauge invariant variable X{r, t), 
ea. (|53|l . Indeed, this equation includes one dimensional delta functions coming from expanding the step function 
in terms of Aro(i) and Ar*(t). Concretely, this tells us that the radial derivative of Gi is not smooth across the 
boundaries, but rather obeys 



lim 



Gi''{ro + e,t)~Ci'iro^e,t) 



lim 



C'/ (r, +e,t)-Gi"^ (r,-e,t) 



-87rG6crAro(t) 
SnGeaAr^t) 



(A29) 
(A30) 



which we will regard as solving for Aro(i) and Ar*(i) once all other functions have been set by the other boundary 
conditions. It is interesting to note that the appearance of one dimensional delta functions in our solutions is in line 
with the arguments recently presented in J5q| . 

Once all boundary conditions have been imposed, we have explicit expressions for all the integration "constants" 
J^i{t) and we also have three additional equations which turn out to give one of the Friedmann equations (involving the 
combination do/ao — (ao/ao)^) and two conservation of energy equations (one for the energy density on each brane). 
Consistency between these three expressions allows us to recover the expression for the Hubble rate (do/oo)^ in terms 
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of the energy densities. The expression written in terms of the six dimensional quantities are not very instructive 
however, so we will not write them out. 

Armed with full solutions, we can now work out the effective four dimensional quantities, which as explained in the 
main text, are given by 

/•"To i-Tr/k-4> 

CT^-*) = 27r / e^^'^^'^adr = 2n e'^°^'''>adr (A31) 

Jo Jir/k-c/i-ro 

p^^\t) ^ 2tt e^«Mp(t)rfr + 27r / e'^«('"VCi(r, i)dr + 27rcre^°(''»' Aro(i) (A32) 

Jo Jo 

fro rro 

p^'^\t) = 2tt e^oMpWrfr- - 27r / e^°WcrCi(r, t)dr - 27r(Te^°('^°)Aro(i) (A33) 

Jo Jo 

p7v/k — (p p-jr/k—ip 

pi'^\t) = 2tt 6^°^''^ p,{t)dr + 2tt e^''WcrCi(r,i)dr 

J TT /k—<p—rQ J Tz/k — cp—rQ 

- 27rCTe^°(''*)Ar,(i) (A34) 

n7T/k — <f) p7r/k — (p 

pi^\t) = 2ti e^«Mp,(t)dr-27r / e^°(''VC7i(r, i)dr 

J TT /k—(p—rQ J 7r/k—(p--rQ 

+ 27rCTe^°(''*^Ar,(i). (A35) 

Plugging these into the expressions for the Friedmann equations and conservation of energy allows us to derive 
eas. (|()7l70t . 
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